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ABSTRACT 

We investigate the theoretical expectations for detections of supermassive binary black holes that can be 
identified as sub-parsec luminous quasars. To-date, only two candidates have been selected in a sample 
comprising 17,500 sourc es selected from the Sloan Digital Sky Survey (SDSS) Quasar Catalog at z < 0.70 
(Boroson & Lauer 2009). In this Letter, we use models of assembly and growth of supermassive black holes 
(SMBHs) in hierarchical cosmologies to study the statistics and observability of binary quasars at sub-parsec 
separations. Our goal is twofold: (1) test if such a scarce number of binaries is consistent with theoretical 
prediction of SMBH merger rates, and (2) provide additional predictions at higher redshifts, and at lower flux 
levels. W e determine the cumulativ e number of expected binaries in a complete, volume limited sample. Mo- 
tivated bv lBoroson & Lauerl d2009). we apply the SDSS Quasar luminosity cut (Mj < —22) to our theoretical 
sample, deriving an upper limit to the observable binary fraction. We find that sub-parsec quasar binaries are 
intrinsically rare. Our best models predict ~ 0.01 deg~ 2 sub-parsec binary quasars with separations below 
~ 10 4 Schwarzschild radii (u or b > 2000 km s -1 ) at z < 0.7, which represent a fraction ~ 6 x 10 -4 of un- 
absorbed quasars in our theoretical sample. In a complete sample of ^10,000 sources, we therefore predict 
an upper limit of ~10 sub-parsec binary quasars. The number of binaries increases rapidly with increasing 
redshift. The decreasing lifetime with SMBH binary mass suggests that lowering the luminosity threshold does 
not lead to a significant increase in the number of detectable sub-parsec binary quasars. 
Subject headings: black hole physics — cosmology: theory — galaxies: nuclei — quasars: general 



1. INTRODUCTION 



SMB Hs appear to inhabit most galaxy centers (e.g., 



iRichstone et al.lll998blFerrarese & Fordl2005l) . and in ACDM 
cosmologies galaxies experience multiple mergers during 
their cosmic assembly. SMBH binaries (SMBHBs) are there- 
fore expected to be recurrent, albeit transient features of most 
galactic bulges. Observationa lly, the paucity of q uasar pairs 
(~0.1%) on galactic scales (Fore man et al.l 12009) points to- 
ward rapid inspiral of the SMBHs down to (sub-)p arsec scales 
where they form a gravi tationally bound pai r (IMayer et al.l 



2007). As discussed in iForeman et alJ (120091) if we assume 
that every galaxy hosts a SMBH and that quasar activity is 
triggered by galaxy mergers, the probability of observing a 
double quasar scales with the ratio of the quasar lifetime to 
the total merger timescale as ~ (i gS oAmerge) 2 , if the two 
quasars light up at different, random times. This is con- 
sistent with the fract ion of quasar pairs found in the SDSS 
(lHennawietaDl2006l), ift qsn < 10 V an d t merge ~ 10 9 yr 
(for mo re details. see | Foreman et al.ll2009l) . At lower level of 
activity, IComerford et all d2008l) find that about 2% of early- 
type galaxies host candidate AGN pairs in the same galaxy. 

Detecting sub-parsec binaries by imaging techniques is 
extremely difficult. However, the presence of SMBHBs in 
AGNs can be discovered spec troscop i cally, a s double broad- 
line emission systems (see iGaskelll ([T996) and references 
therein). The two sets of broad emission lines originate in 
gas associated with the two SMBHs, and the velocity sepa- 
ration between the two emission line s ystems tra c es the pro- 
jected orbital velocity of the binary. iDotti et ail (120081) and 
Bogdanovi c et al.l ((2008) extended the pre-existent spectro- 
scopical technique, discovering that, depending on the SMBH 
mass ratio, the AGN spectrum shows two sets of broad lines 
(equal mass binaries) or a single set of broad lines and two 
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sets of narrow emission lines at different redshifts (unequal 
mass binaries). In the latter case, only one of the two SMBHs 
is active, and the two sets of narrow emission lines correspond 
to emission from low density gas in the potential well of the 
binary and from the "standard" narrow line region (NLR) of 
th e AGN. 

IDotti et al l (120081) and iBogdanovic et all (120081) ap- 
ply the binary model to the peculiar quasar SDSS 
J092712.65+294344.0. This quasar exhibits two distinct sets 
of lines. The first set of very narrow emission lines is assumed 
to be emitted in the NLR of the host and trace the redshift of 
the host galaxy (z = 0.713). The second set comprises two 
blue-shifted systems featuring different FWHM: the broad 
Mg II and Balmer emission lines with FWHM » 4000 km 
s -1 , and narrow lines with FWHM « 460 - 2000 km s _1 , 
both consistent with a redshift of 0.698. In this model, the 
source emitting the blue-shifted line system is gas inside or 
near the broad line region (BLR) of the secondary, comoving 
with the SMBH with a light-of-sight velocity of 2630 km s _1 
re lative to the rest-frame o f the host. 

Bor oson & Lauerl d2009) developed a principal components 
analysis technique that identifies sources having peculiar 
spectral characteristics. They applied this procedure to the 
restframe optical spectrum of ^17,500 quasars. Their sam- 
ple comprises all quasars having z < 0.7 from the fifth 
release of the SDSS Quasar Catalog (ISchneider et al.ll2007l) 
plus all sources classified as quasars in the seventh SDSS 
data release. Of the 17,500 objects in their entire sam- 
ple, only two objects have multiple redshift systems con- 
sistent with the presence of a candidate SMBHB, SDSS 
J0927 12.65+294344.0, a nd SDSS J153636.22+044127.0 (but 
see Chornok et al. 2009; iHeckman et"all l 2009l: IShields et all 



t et alJ|20 
l200§~ 



2008; IWrobel & Laodl2009t iGaskell 

In this paper we assess the expected number of merging 
SMBHBs based on realistic merger rates of SMBHs in hier- 
archical cosmologies, and determine an upper limit to sub- 
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parsec binary quasars detectable as double broad-line emis- 
sion quasars 

2. BLACK HOLE MERGER RATE AND QUASAR ACTIVITY 

We trace the evolution of SMBHs within a plausible sce- 
nario for the hierarchical assembly, growth, and dynamics 
of SMBHs in a ACDM cosmology. Our model has been 
shown to capture many features of the SMBH population 
(e.g., luminosity function of quasars, X-ray background, 
SMBH mass density). The m ain features of the models 
have been discussed el s ewhere dVolonteri et ail I2003L 120051 : 
IVolonteri &Reesl 120061: ISesana et al.l (2007b. and references 
therein). We summarize in the following the relevant assump- 
tions. 

SMBH "seeds" form at high redshift (z > 15) in highly 
biased halos, co rresponding to 3.5-4 a peaks of the density 
fluctuation field (IVolonteri et alJ2007l) . The initial SMBH oc- 
cupation fraction is therefore low, but along the cosmic hi- 
erarchy SMBHs are incorporated into massive and massive 
systems, as galaxies grow in a ACDM cosmology. 

The occupation fraction of SMBHs increases with time, 
and approaches uni t y for massive galaxi es at low redshift 
(iMarulli et al.l 120061: IVolonteri et all 120071) . In our scheme, 
therefore, not all galaxy mergers lead to SMBH merg- 
ers, but only those involving two galaxies both hosting 
SMBHs. We further assume that SMBHs merge within the 
merger timescale of their hosts, which is a likely assump- 
tion fOTjSMTiHbin^es^^ merg- 
ers dUcaEEiDllQOl HI! |Diditan|200l,l003). This 
is the most likely scenario in the context of this work, as 
quasar fueling requires a substantial gas supply. We explored 
an alternative scenario where at late cosmic times SMBHBs 
shrink via three-body interactions, i.e., by capturing and eject- 
ing at much higher velocities the stars passing by within a 
distance comparable to the binary separa tion dMerritll 120061: 
IVolonteri et al.|[2003t ISesana et al.ll2007al) . and we found that 
the merger rate in the z < 1 redshift range is very simi - 
lar, as already fou nd in previous tests (ISesana et al.l 12005b . 
SesanaetaL|(2008) compare our theoretical merger rates with 
the merger rates inferred by observations of the fraction of 
close galaxy pairs (assuming that SMBH masses scale with 
bulge masses). Our SMBH merger rate is consistent within 
a factor of < 2 with the merger rate of massive spheroidal 
found in Bell et al. (2006), who quote a factor of 2 uncer- 
tainty in their rate estimate. 

We base our model for SMBH mass growth on a set of sim- 
ple assumptions, sup ported by both simu lations of AGN trig- 
gering and feedback ( Springel et al. 2005), and analysis of the 
relationship betwe en SMBH masses (Mbh) and the proper- 
ties of their hosts dGebhardtet all 120001: iFerrarese & Merrittl 
2000; iFerraresel 120021) . SMBHs in galaxies undergoing a 
major merger (i.e. having a mass ratio > 1 : 10) un- 
dergo accretion. Each SMBH accretes an amount of mass, 
AAf = 9 x 10 7 M Q (cr/200kms~ 1 ) 4 , that scales with the 
A^bh — c* relation of its hosts (see Volonteri & Natarajan 
2009). Accretion starts after a dynamical timescale and lasts 
until the SMBH has accreted AM. 

The accretion rate during the active phase is derived 
from the empirical distr ibution of Eddington ra tios, A = 
log(Lboi/iEdd), found in Me rloni & Heinzl(l2008l) . We adopt 
a fitting function of the Eddington ratio distribution as a func- 
tion of SMBH mass and redshift (Merloni 2009). The dis- 
tributions of Eddington ratios, f\, are computed in 10 red- 
shift intervals (from z = to z = 5) for 4 different mass 



bins (6 < log(M BH /M Q ) < 7, 7 < log(M BH /M Q ) < 8, 
8 < log(M BH /M Q ) < 9, 9 < log(M BH /M ) < 10), 
and then fit with an analytic function which is the sum of a 
Schechter function and a log-normal. The Eddington ratio 
distributions are normalized so that at a given mass and red- 
shift, J x f\ = 1. We adopt the bolometric corrections pre- 
sented in [Richards et al.l (120061) , and we correct for ab sorbed 
quasars according to model 4 in lLa Franca et al.l {2005). 

Our first check is on the cumulative number counts of 
quasars, regardless of binarity. We select all accreting SMBHs 
in our theoretical sample, and apply the same luminosity cut 
as in the SDSS Quasar Catalog, M 4 < -22. This choice 
is motivated by the SDSS being the largest quasar catalog 
currently available, however we stress here that our models 
provide upper limits to the number of observable binaries, in 
a complete volume limited sample. We compare our num- 
ber counts to the expected number from the integration of 
the bo lometric luminosity function of quasars (Hopki ns et al.l 
2007) at Mi < -22. The comparison is shown in the top 
panels of Fig. 1 (top histogram and dashed curve). 

3. SUB-PARSEC BINARY QUASARS 

To evaluate the expected number of binary quasars, we start 
from the SMBHB merger rate. Our theoretical models indi- 
cate which of these binaries are in an active phase, following 
a major merger. We initially analyze the complete sample of 
SMBHBs, that is, we consider all merging SMBHs, regardless 
of their activity status. This first model (model I) provides a 
strong upper limit to the number of theoretical sub-parsec bi- 
nary quasars. To each SMBH in our sample we randomly 
assign an Eddington ratio, A, from the normalized distribu- 
tion (Merloni 2009). Note that, since -5 < A mm < -3 for 
most masses/redshift this model does not allow for 'quiesc ent' 
SMBHs at, e.g., the level of Sgr A* dOuataert et al.ll 19991) . 

We further assume that only SMBHBs with a mass ratio, 
q = M B h,2/M B h,i ^ 1, above a certain threshold (q > 10 _1 
or q > 10~ 2 ) create distinguishable double broad emission 
line systems. We regard q > 10 _1 as our best choice, be- 
cause of an addition al independent motivation, as follows. 
Calleg ari et~aT1 (f2008) study the formation of SMBHBs during 
galaxy mergers, and how SMBH pairing depends on the in- 
terplay between different physical processes (dynamical fric- 
tion, tidal and ram-p ressure stripping). From this analysis, 
Calle gari etaTI (120081) find that q = 10 _1 is the minimum 
mass ratio between two merging galaxies (and as a conse- 
quence between the two SMBHs) to guarantee the formation 
of a SMBHB. However, tidal and gas-pressure stripping can 
be reduced for extreme structural properties and orbital pa- 
rameters of merging galaxies. We consider the very conser- 
vative case of SMBHBs with q > 10 -2 to take into account 
every possible merger configuration (e.g., plunging radial or- 
bits). 

Finally, we assign a lifetime to SMBH binaries detectable 
at sub-parsec separations following the detailed study of 
SMBHB d ynamical evolution in circumbinary disks per- 
formed by Hai man et alj (120091) . The lifetime corresponds to 
the time spent by the SMBHB at a separation such that the 
shift between the BLRs (or the NLR and the BLR) is a few 
thousand kms^ 1 , b ecoming compar able with the widths of 
typical broad lines (She n et aTll2008l) . Blending of profiles 
with smaller velocity differences would be missed (Gaskell 
1996). For v or b > 2000 km s" 1 , the typical separation is 
r = 11.25 x 10 3 Schwarzschild radii, and for these condi- 
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tions, lHaimanet all ([2009) give a lifetime: 



tlife = 6Myr 
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where Mbm = Mbh,2 + Mbh,i, and we have chosen the 
longest residence time, thus providing an upper limit to SMB- 
HBs lifetime. 
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do produce distinguishable double broad emission line sys- 
tems, then the expected SMBHB merger rate must be lower 
than we predict. The high luminosity selection criterion leads 
to a sample composed of actively accreting (A ~ —1) massive 
(~ 10 8 M ) SMBHs. 

In a second model (model II) we select only SMBHBs 
where at least one of the SMBHs is active, according to our 
merger-driven quasar activity scheme. The results of this 
model are shown in the top panels of Fig. 1. In this case 
the expected number of binaries per deg 2 is 0.02+QQ2 f° r 
q > 10" 2 and O.Oltg;^ for q > 10 _1 . Detectable binaries 
represent a fraction < 6 x 10~ 4 — 10~ 3 of the unabsorbed 
quasars, consistent with Boro son & Lauerl d2009) findings (a 
fraction ~ 1 - 2 x 10" 4 ). 

If we decrease the luminosity threshold we expect two 
factors enter into play: on the one hand, the merger rate 
of SMBHs is expected to increase at lower SMBH masses , 
where the mass function is less steep (Gul tekin et al.ll2009l) . 
This is indeed what we find in our models (Fig. 12 see also 
Dotti et al. 2009 and Sesana et al. 2005). On the other 
hand, however, the lifeti me of detectable b inaries decreases 
with decreasing mass (Haim arTet al.l 12009). making the de- 
tection harder. Using the scaling for lifetimes presented in 
Haiman et al. 2009, we indeed expect that the number of 
detectable sub-parsec binary quasars does not increase dra- 
matically with decreasing luminosity, because of the shorter 
timescale over which they are observable. For instance, if we 
decrease the flux limit by a factor of ten, we find negligible 
changes in model II, and a mild increase in the number of bi- 
naries by a factor of 3 in model (I): 0.17 ± 0.13 for q > 10 _1 
and 0.79 ± 0.64 for q > 10~ 2 . 



FIG. 1. — Cumulative number of sub-parsec binary quasars with Mi < 
—22, per square degree. Left panels: q > 10~ 2 . Right panels: q > 10 _1 . 
Bottom: we assume that each SMBH is at some level "active", i.e. we do 
not consider a merger-driven quasar activity (model I). Top: we select only 
SMBHBs in galaxies that have recently experienced a major merger, i.e. we 
do assume that quasar activity is merger driven (model II). In the upper-left 
panel we compare our total quasar sample (top curve with error-bars) to the 
expected number co unts from the integr ation of the bolometric luminosity 
function of quasars ( Hopkins et al. 2007) at M, < -22 (dashed curve). In 
the upper-right panel we also show the number of unabsorbed quasars when 
we correct for absorption (top curve with error-bars). 

The results of this model (model I) are shown in Fig.[T](bot- 
tom panels). We find that at z < 0.7 the total expected num- 
ber of binaries with Mi < -22 per deg 2 is 0.27 ± 0.27 for 
q > 10~ 2 and 0.041°;°^ for q > 10" 1 . These detectable 
binaries represent a fraction ~ 10~ 2 and ~ 2 x 10~ 3 re- 
spectively of the unabsorbed quasars. Within t he sam ple of 
17,500 sources analyzed by IBoroson & Lauerl 12009), 9895 
objects, including the two putative binary SMBHs, belong 
to the uniformly selected statistical sample of SDSS quasars 
(Boroson, private communication). The statistical sample 
reaches a completeness larger than 90% at z < 1 for sources 
with apparent magnitude i < 19.1, roughly corresponding to 
Mi < 24 at z < 0.7 (Richards et al. 2002). Given these 
caveats, our results must be considered upper limits to the 
number of detectable sub-parsec binary quasars. We there- 
fore find margina l agree ment, within the uncertainties, with 
IBoroson & Lauerl (120091) findings 2 . If binaries with q < 10 ~ 2 

2 We have calculated the error by using the simple Poisson statistics, 1— cr 
confidence. Th e errors are lower limits given the small number statistics. See 
IGehreldfi986l) . 



4. DISCUSSION 

Stimulated by the recent putative discovery of two 
candidate sub-parsec SMBHBs identified as quasars with 
multiple redshift line systems (BLR and NRL), SDSS 
J0927 12.65+294344.0 and SDSS J153636.22+044127.0, we 
investigate theoretically the occurrence of sub-parsec SMB- 
HBs that can be identified as sub-parsec binary quasars. We 
study the SMBH cosmic evolution via a Monte-Carlo merger 
tree approach. We trace the growth and dynamical history 
of SMBHs from high redshift via physically motivated pre- 
scriptions, that allow us to reproduce many observational con- 
straints. 

Our approach provides us with a catalog of SMBHBs, 
for which we know the masses and the redshift. We fur- 
ther assume that only SMBHBs with a mass ratio, q = 
Mbr,2/Mbh,i ^ 1, above a certain threshold (q > 10 _1 or 
q > 10~ 2 ) create distinguishable double broad emission line 
systems. Motivated by the work by Boroson & Lauer (2009), 
we apply to our theoretical sample the same luminosity cut 
as in the SDSS Quasar Catalog, Mi < —22, thus deriving 
upper-limits to the fraction of detectable binaries. We stress 
here that the theoretical sample is complete and volume lim- 
ited. Since the S DSS quasar catalog is not a complete sample 
dSchneider et al.ll2"007l) . a direct comparison with Boroson & 
Lauer (2009) is not appropriate. However, our results are con- 
sistent with the binary fraction derived for the subset of the 
17,500 quasars used by Boroson & Lauer that are part of the 
statistical sample of the SDSS. 

Our merger-driven quasar scheme provides us also with an 
accretion rate, hence a luminosity. We analyze two models 
that likely bracket the theoretical uncertainties. Model I ig- 
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FIG. 2. — Merger rate of SMBHs per unit redshift. Top histogram: binary 
mass > 3 X 10 7 Mq. Bottom histogram: binary mass > 3 X 10 8 Mq. 

nores merger-driven quasar activity and we assume that each 
SMBH is at some level "active". Each SMBH in our binary 
sample is assigned an Eddington ratio, A, from the normalized 
distribution derived from synthesis model for AGN evolution 
dMerloni & Heinzll2008l) . Model I is therefore our strong up- 
per limit to the number of detectable SMBHBs. Model II 
is more rooted into our quasar activity scheme, as we select 
SMBHBs where at least one of the SMBHs is active, accord- 
ing to our merger-driven scenario. 
Our main findings are as follows: 

• Sub-parsec binary quasars are intrinsically rare, due 
to a combination of strict requirements: the time over 
which SMBHB are detectable through line shifts de- 
creases with decreasing binary mass. On the other 
hand, the merger rate of SMBHBs increases with in- 
creasing mass. 

• In a volume limited, complete sample of ~ 10,000 
sources at z < 0.7, our best models (II), that re- 
late quasar activity to galaxy mergers, predict an up- 
per limit of ~ 5-10 sub-parsec binary quasars. Model 



I, that does not associate qu asars to mergers, i s only 
marginally compatible with Boroson & Lauerl |2009) 
who find only 2 candidate sub-parsec binary quasars in 
the statistical SDSS quasar sample. 



• Fig. 1 extends our predictions out to z = 1. The number 
of detectable binaries increases by a factor ~ 5 — 10 
from z = 0.7 to z = 1. 



• The lifetime over which SMBHBs can be detected as 
sub-p arsec quasars decre ases with decreasing binary 
mass tHa iman et al.ll2009l) . This effect is stronger than 
the increase in the merger rate of SMBHB at lower 
masses. Lowering the luminosity threshold is unlikely 
to lead to a large increase in the number of detectable 
sub-parsec binary quasars. 
SDSS-III will increase the spectroscopic quasar sample 
and will provide a good testing ground for our predictions. 
The masses of SMBHBs that can be identified as sub-parsec 
binary quasars are too large for the gravitational waves 
emitted by these binaries to be detectable at merger by the 
Laser Interferometer Space Antenna (LISA), which will 
instead focus on the mass range 10 5 — 1O 7 M . However, 
such massive binaries (in a later evolutionary stage, when 
the binary has shrunk by an additional factor of ten and the 
dynamical evolution is driven by emission of gravitational 
radiation) are typical candidates for detection via Pulsar 
Timing Arrays (PTAs, e.g. the Parkes radio-telescope). PTAs 
rely on the effect of gravitational waves on the propagation 
of radio signals from a pulsar to the Earth, producing a 
characteristic signat ure in the time of arrival of radio pulses. 
ISesana et al.1 (|2008) find that the mass distribution of the 
SMBHBs detectable via PTAs peaks at — 10 8 M Q , with most 
binaries at z < 1, the same mass range probed by sub-parsec 
binary quasars identifiable in the SDSS. 
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